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G\ \ Abstract 



Ground state of the bond-random antiferromagnetic S = 1 Heisenberg chain with the biquadratic 
Q ! interaction — PJ2t(^i ' Si+i) 2 is investigated by means of the exact-diagonalization method and 
qq | the finite-size-scaling analysis. It is shown that the Haldane phase (3 ~ persists against the ran- 
^NJ domness; namely, no randomness-driven phase transition is observed until at a point of extremely- 
broad-bond distribution. We found that in the Haldane phase, the magnetic correlation length is 
\ kept hardly changed. These results are contrastive to those of an analytic theory which predicts 
' a second-order phase transition between the Haldane and the random- singlet phases at a certain 
^ critical randomness. 
O > 

1 Introduction 

The ground-state phase transition of the disordered quantum system has attracted much attention 
recently Because the transition is driven solely by the quantum fluctuation, the transition 
r lies in different universality from the finite-temperature transitions. The randomness contributes 
significantly to the anisotropy of the real-space and the imaginary-time directions. Thereby, the 



O randomness would also promote new types of phase transition. 

For this purpose, the ground state of the random quantum spin system has been studied 
^ ■ considerably, see the article fl], for a review. In one dimension, in particular, a real-space- 
^ decimation procedure || f|] yields various definite predictions |@, 0, 0|. For the Heisenberg model 
in the presence of the bond randomness f| @, H, for example, the theory tells that infinitesimal 
randomness drives the ground state to the random-singlet phase, where the averaged ground-state 
correlation decays obeying the power law, 

[<^>L ~ Vr 2 • (i) 

Note that the correlation decays faster than that of the pure chain ~ 1/r. Numerical simulations 



followed [|K], [11], |TI| so as to confirm the analytic predictions. It would be noteworthy that the 
5=1/2 spin system is viewed as a hard-core boson system ||13|| . The random spin system is, thus, 
equivalent to the strongly repulsive boson system embedded in a certain randomness, which is of 
the current interest |T3| . 

In these studies, the magnitude of the constituent spin is set to be one half. In one dimension, 
however, there exist some other sorts of magnets such as the integer-spin Heisenberg chain and 
the 5 = 1/2 Heisenberg ladder model. The ground states of these magnets are different from 



the ground state of the 5 = 1/2 Heisenberg chain |TJ|, [TH] : In the former case, a finite magnetic 



excitation gap opens above the ground state, and the magnetic correlation decays exponentially in 
the ground state. In the present paper, we investigate the bond-random antiferromagnetic 5 = 1 
spin chain. The Hamiltonian is given by 

n = J i {Si • s !+ i - p(Si ■ s m ) 2 } . (2) 



The operators Sj = (Sf, Sf, Sf) denote the S = 1 spin operators acting on the site i, and satisfy 
the periodic boundary condition, S^+i = Si. We used the following bond distribution, see Fig. |T], 
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Figure 1: The probability distribution of the random bond J. The randomness ranges as < A < 1. At A = 1, 
the distribution becomes extremely broad; namely, infinitesimally weak bond appears. 

P A (J) = i(6(J-l+A)-e(J-l-A)) . (3) 

(The function Q(x) denotes the step function.) Note that only the antiferromagnetic bond appears. 

Because of the presence of the excitation gap, the system (0) lies beyond the scope of the 
previous studies of the S = 1/2 magnets explained above. It is quite suggestive that the real- 
space decimation theory becomes inapplicable for the cases other than 5 = 1/2. In order to adapt 
the decimation procedure even for the case S = 1, several authors proposed schemes to map the 
random S — 1 chain to an effective 5 = 1/2 model [ I7[ , |I8 |. Their theory and the consequences 
are reviewed afterwards. 

The effect of the non-magnetic impurity doping upon the spin-gap ground state has been 



considered rather throughly |L9], |2(J, ^TJ: The non- magnetic impurity causes notable influences for 
the Heisenberg ladder. Infinitesimal doping distracts the spin gap, and gives rise to a magnetically 
(quasi-) ordered ground state with gapless excitation. The effect of the bond randomness, which 
is of the present interest, would be more subtle, and remained unsolved. 

For the pure system { Jj = 1}, the ground-state phase diagram of the Hamiltonian (^) is known, 



see the paper [27 for details. We explain the region (3 ~ relevant to the present study. In the 
region —1</3<1, the system is in the Haldane phase, which is mentioned above; namely, a finite 
excitation gap opens, and the magnetic correlation decays exponentially in the ground state. At 
the critical point (3 = 1, the Hamiltonian is integrable. It is predicted in a field-theoretical manner 
[P3"| that the criticality is of the central charge c = 3/2, and the biquadratic interaction generates 
the excitation gap in the form, 

AE oc IT , (4) 

where v — 1. Numerical simulations have tried to confirm the above scenario E7J. 
A strong finite-size corrections, however, arise so as to prevent definite conclusions. Typical 
simulation claims that a transition locates at (5 « 0.5, and gapless phase extends over 0.5 < (3 



2R EH]. The inconsistency might be originated in the long correlation length £ in the Barber- 



Batchelor phase 1 < /?; £ ~ 21 at ft — ► oo p8| . The length apparently exceeds the limit tractable 
by means of the exact-diagonalization method. 

The paper is organized as follows. In the next section, we explain the previous analytic theory 
employing the real-space-decimation method ||T7| , [I~8|j . The theory predicts that a phase transition 
separates the Haldane and the random-singlet phases at a certain critical randomness. In the 
section ^, we present our result by means of the exact diagonalization method. We show that the 
finite-size scaling assuming the theoretical prediction fails; namely, the present numerical estimate 
of the transition point and the critical exponent contradicts the prediction. The Haldane phase 



persists against the randomness until the extremely-broad-bond distribution A = 1. The results 
of the correlation length and the spin stiffness confirm the present conclusion. In the last section, 
we summarize our results, and discuss a possible scenario interpreting our numerical results. 



2 Review of the predictions with the real-space decimation analysis 



In this section, we review the predictions of the analytic theory |17|, |18fl . Some are contrastive to 
our numerical results shown later. As was introduced in the section [3], the real-space-decimation 
method describes the one-dimensional random magnets quite successfully. The method, however, 
fails in the cases other than S = 1/2. Hyman and Yang thus proposed the mapping with which 
the S = 1 random chain (^) is transformed to an effective random S = 1/2 chain. The resultant 
S = 1/2 chain constitutes the alternating two types of random bonds, {J21} and {J2i-i}- The 
even bonds are either ferromagnetic or antiferromagnetic, whereas the odd bonds are always 
antiferromagnetic. 

We sketch how the effective 5 = 1/2 random chain is derived. Suppose a sector within the 
S = 1 chain in which sector the deviation of the randomness is incidentally small. They assumed 
that such sector could be replaced with that without any randomness. As a consequence, the 
random S = 1 chain consists of finite uniform S = 1 sectors coupling antiferromagnetically. 
On the other hand, it was reported that at an end of the uniform S = 1 open chain, the S = 
1/2 magnetization appears spontaneously with a certain localization length [29]. Both the end 



magnetizations couple antiferromagnetically (ferromagnetically) , if the length of the open chain 
is even (odd). Thereby, the finite uniform 5* = 1 sector is replaced with the two S = 1/2 spins 
coupling either antiferromagnetically or ferromagnetically. 

The real-space decimation procedure is employed afterwards. They found a fixed point at a 



critical randomness A c , which separates the Haldane phase and the random- singlet phase [ ^q| . 
(Refer to the section [l] for these phases.) Furthermore, the theory concludes that the correlation 
length diverges as 

£~1/|A-A c r (5) 

with v « 2.3. 

The decimation procedure was simulated numerically as well 0]. The simulation shows that 
an intermediate phase appears between the Haldane and the random-singlet phases. 



3 Numerical results 



In this section, we investigate the Hamiltonian (|2|) by means of the exact-diagonalization method. 
The exact-diagonalization method has been utilized successfully in the course of the studies of the 
1/2 random spin systems 



S 



3.1 Suppression of the string correlation by the randomness 



In Fig. 0, we plotted the string correlation |32], |33| 0* tT - (L/2) against the randomness A at the 
Heisenberg point (3 = 0. The string correlation is defined as 



^string (i V 



(£*e i7r £t-^) 



(6) 



The random average [ ] av is taken over ninety samples throughout the present study. The string 
correlation was invented in order to detect characteristics of the Haldane state. The long-range 



limit was estimated numerically, 0" t 



string V 



00) -> 0.38, at (3 = 



(It is noteworthy that the 



usual correlations, such as the Neel correlation, are short-range.) In Fig. |2|, we see that the string 
order dominates in the weak-random region, while it becomes suppressed as the randomness is 
strengthened. 




Figure 2: The string correlation 0^ tlins (L/2) (|^) is plotted for various randomness A at the Heisenberg point 
(3 = 0. The string correlation which characterizes the Haldane phase becomes suppressed as the randomness is 
strengthened. 



In order to observe whether the string correlation is long-range or not, we evaluate the Binder 



parameter [35 for the string correlation, 



where O = Z)j=i e 17r ^ fc = 1 k S*. The Binder parameter is enhanced (suppressed) with increasing 
the system size, if the corresponding order is long-range (short-range). It is plotted in Fig. [| 
where the parameter varies in the same range as that in Fig. |2|. We see that the string order 
develops in the whole range of the randomness. The result contradicts the picture |TJj reviewed 
in the previous section. 

As the biquadratic term is turned on (/3 > 0), the Haldane phase becomes suppressed. (Note 
the phase diagram for the pure case reviewed in the section [1].) In Fig. f|, where we plotted the 
Binder parameter at (3 = 0.2. It is implied that there exists a transition point (3 ~ 0.6; in one 
side (3 < 0.6 the string order develops, while in the other side 0.6 < (3 it is disturbed. Some might 
wonder that the result of (3 = 0.2 in Fig. [| is inconsistent with the result of (3 = in Fig. [|, 
and the former supports the picture of the previous theory [[Ti|. In the following, however, we 
show an evidence that the transition point (intersection point) is extrapolated to A m 1 in the 
thermodynamic limit L — > oo. 

3.2 Analysis of the transition point 



According to the theory [18] reviewed in the section 0, the transition should be the second order. 
Namely, the Binder parameter ([7]) obeys the finite-size scaling hypothesis, 

U(L) = L X U ((A - A c )L 1/,y " 



The scaling dimension x is zero, because the Binder parameter is dimensionless. Therefore, the 
data (A — /\ c )L l ^)-U (L), the so-called scaling data, must align along a curve which is independent 
on the system size L. We determine the scaling parameters A c and v in eq. (|5p, so that the data 
converge to such an universal curve; see the appendix for details. 

In Fig. |5|, we show the scaling plots for (3 = 0.2. In order to see the systematic corrections to 




Figure 3: The Binder parameter U(L) ([?]) is plotted for the same parameter range as that shown in Fig. ||. 
The plot shows that the string order develops over the whole region, and at the extreme randomness A = 1, it is 
disturbed. 




Figure 4: The Binder parameter U(L) is plotted for various randomness A at (3 = 0.2. We see an intersection 
point A ~ 0.6 which suggests a transition point. 
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Figure 5: Scaling plots for /3 = 0.2. In order to observe the systematic corrections to the finite size scaling, we 
managed to analyse different three sets of the system sizes; (a) L = 6, 8, 10, (b) L = 8, 10, 12 and (c) L — 10, 12, 14, 
respectively. 



the finite size scaling, the scaling analyses are managed for different three sets of the system sizes; 
(a) L — 6, 8, 10, (b) L — 8, 10, 12 and (c) L = 10, 12, 14, respectively. We define the approximative 
transition point A(Z) and the exponent which are determined by the scaling plot L = I — 2, / 
and 1 + 2. We found that the scaling parameters result in the following, 

A c (8) = 0.67, A c (10) = 0.70, A c (12) = 0.87, 
i/(8) = 0.92, i/(10) = 0.76, i/(12) = 0.84. (9) 

The transition point converges to A ~ 1 for large system sizes, and the exponent approaches 
v ~ 0.9. It is suggested that the transition point is identical with that at the Heisenberg point 
(3 = 0, see Fig. g 

We show the scaling plots for f3 = 0.3 as well, see Fig. §. The plots yield the following estimates: 

A c (8) = 0.33, A c (10) = 0.47, A c (12) = 0.56, 
i/(8) = 0.65, i/(10) = 0.94, i/(12) = 0.34. (10) 

We observe that the transition point approaches A — > 1, while the exponent rather scatters. The 
scaling analyses beyond (3 ~ 0.3 are suffered by strong finite-size corrections which also exist 
for the uniform system (refer to the section [l]). In consequence of the previous and the present 
subsections, we observe that a phase transition might occur at A = 1. The universality class 
is rather ambiguous. It is suggested, at least, that the previous prediction v w 2.3 Jl8| is not 
accordant with our numerical analysis. 

In the next subsection, we show an evidence that the transition is the first order. Even for 
such case, the above finite-size-scaling analysis is still meaningful in the sense that it yields the 
estimate of the intersection point of the Binder parameters. The intersection point is expected to 
converge to the transition point even for the first-order transition. The exponent, on the contrary, 
is not universal, and possibly would not be well-defined for very large system sizes. 

3.3 Correlation length and the spin-stiffness constant 

In order to confirm that a first-order transition takes place at A = 1, we show the result of the 
correlation length. We found that the correlation length is remained unchanged over the whole 
region < A < 1; namely, the correlation length does not diverge as in eq. (|5|). 
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Figure 6: Scaling plots for = 0.3. The analysis is the same as that shown in Figs. [|. 

In Fig. |7|, the averaged Neel correlation is plotted for various randomness at the Heisenberg 
point (3 = 0. The averaged Neel correlation is given by 



0Neel(j'-*)= <(-ir^> 



The decay rate ~ q-U- 1 )/ 6 - 2 estimated previously [^] at A = is shown as well. The magnetic 
correlation length is hardly changed with the randomness. Moreover, we observe that the corre- 
lation itself is not changed very much. The result indicates that the Haldane phase with finite 
correlation length £ ~ 6.2 continues up to the extreme point A = 1. 

In Fig. ||, the spin stiffness constant is plotted. The stiffness constant |37| is defined by 



Ps 
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where 9 denotes the boundary twist, 
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The stiffness constant 
ness destroys long-ran^ 
vanishes exponentially 
magnetic correlation, 
the while randomness 
in the Haldane phase, 
it transits at the point 



is generally employed for investigating the random system: If the random- 
^e coherence, the stiffness should vanish. In the Haldane phase, the stiffness 
with enlarging the system size, because the ground state has short-range 
We observe in Fig. |8] that the stiffness actually decays exponentially over 
with the localization length unchanged. In consequence, we conclude that 
the ground state property is kept unchanged against the randomness, and 
A = 1 abruptly. We will discuss the transition scenario in the next section. 



4 Summary and discussions 

We summarize and discuss our numerical results. The S = 1 Heisenberg chain (0) with the bond 
randomness A is studied numerically. The Haldane phase (/? ~ 0) survives over the whole random 
region < A < 1: The intersection point of the curves of the Binder parameter converges to the 
point A w 1 in the thermodynamic limit. In fact, the result is consistent with the previous theory 
[p8| , which suggests that the Haldane phase would be stable against randomnesses considerably, 



because the Haldane state has such a hidden correlation as the string correlation (pi). (Hence, the 




Figure 7: Logarithm of the Neel correlation 0^ eol (L/2) is plotted against the system size L for various randomness 
at the Heisenberg point /3 = 0. The decay rate (~ e _r / 6 2 ) at A = Q is shown by the dashed line. The plot 
shows that the correlation decays exponentially in fact, and the correlation length is not changed very much with 
increasing the randomness. 




Figure 8: Logarithm of the spin-stiffness constant p s (L) is plotted against (a) L and (b) logL, respectively, for 
various randomness at f3 = 0. 



theory [TE\ reviewed in the section [2] somehow contradicts the above.) The present system is, 
thus, quite contrastive to the Heisenberg ladder with the non-magnetic-impurities, whose spin- 



liquid ground state is unstable against infinitesimal doping [[19], g0| £1 



We observed that the correlation length is hardly changed with the randomness increased. 
Hence, we conclude that the correlation length does not diverge in the vicinity of A = 1, and 
thus the phase transition is the first order. (Note that according to the theory [18| , a second-order 
transition separates the Haldane phase and the random-singlet phase at a critical randomness A c .) 
At the transition point A = 1, namely, the extremely-broad-bond distribution, infinitesimally weak 
bonds appear. It is thereby expected that at this very point the random-singlet phase is realized 

B 

Finally, we mention what phase extends beyond the point A > 1. In the region, ferromagnetic 
bonds appear as well. This sort of randomness is considered for the 5=1/2 Heisenberg chain, 
see the article jE]] for a review. We suspect that in this region, the ground state is similar to that 
realized for the the S = 1/2 chain: It is found that the low-temperature properties resemble those 
of the classical counterpart [fH| . This similarity is quite convincing, because the ferromagnetism 
is not affected very much by quantum fluctuation. It is natural that at least the magnitude of the 
spin does not alter the ground-state property. The Haldane phase A < 1, on the contrary, owes 
the stabilization to the zero-point quantum fluctuation. In consequence, we see that at A = 1 a 
drastic transition occurs from the Haldane phase to the above rather classical phase. 
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A Details of the present scaling analyses 

We explain the details of our finite-size-scaling analyses, which we managed in the section |3.2| 
in order to estimate the transition point A c and the exponent v. We adjusted these scaling 
parameters so that the scaled data, such as those shown in Figs. |5] and |], form a curve irrespective 
of the system sizes. In order to see qualitatively to what extent these data align, we employ the 
"local lineality function" S(A C , v) defined by Kawashima and Ito |]4"2"fl : Suppose a set of the data 
points {(xi,yi)} with the errorbar {di{= Si/i)}, which we number so that Xj < x i+ i may hold for 
i — 1, 2, ■ ■ ■ , n — 1. For this data set, the local-lineality function is defined as 

n-l 

S = ^2w(xi,y i ,di\xi^ 1 ,y i ^ 1 ,d i - 1 ,x i+ i,yi +1 ,di + i). (14) 

i=2 

The quantity w(xj, y jy dj\xi, yi,di,x k ,y k ,d k ) is given by 



(16) 



A 

where 

— \^k •Ej/Vi y^i ^j)Vk 



and 

A 2 = d) + (^—^d) 2 + (^^d k ) . (17) 

V X k Xi / \ X k X{ / 

In other words, the numerator yj — y denotes the deviation of the point (xj,yj) from the line 
passing two points {xi,yi) and {x k ,y k ), and the denominator A stands for the statistical error of 
(.Hi — y)- And so, w = {{x/i — y)j A) 2 shows a degree to what extent these three points align. 



As the number of the data points n is increased, the statistical error of S would be reduced. 
The corrections to finite-size scaling might increase instead. In the present analyses, we used 
twenty data in the vicinity of the transition point A c . An example of the plot S is shown in Fig. 
H We observe the minimum at v « 0.76, which yields the estimate of the exponent. 




Figure 9: The local linearity function S (|T^) for the scaling data in Fig. |5| (b) is plotted for various v(l — 10) with 
A c = 0.70 fixed. The location of the minimum yields the estimate of the exponent v. 
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